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Abstract 

Accurate determinations of the MS 6-quark mass mi, (mi,) from a (e^e^ hadrons) experimental data cur- 
rently contain three comparable sources of uncertainty; the experimental uncertainty from moments of this cross- 
section, the uncertainty associated with Us (Mz), and the theoretical uncertainty associated with the renormal- 
ization scale. Through resummation of all logarithmic terms explicitly determined in the perturbative series 
by the renormalization-group (RG) equation, it is shown that the renormalization-scale dependence is virtually 
eliminated as a source of theoretical uncertainty in mb (nib). This resummation also reduces the estimated effect 
of higher- loop perturbative contributions, further reducing the theoretical uncertainties in mb {mb). Further- 
more, such resummation techniques improve the agreement between the values of the MS 6-quark mass extracted 
from the various moments of R{s) = a (e"'"e~ — > hadrons) /apt [apt = 47ra^/(3s)], obviating the need to choose 
an optimum moment for determining (mb). Based on this analysis, the resulting value of the 6-mass is 
mb (mb) = 4.207 GeV ±40 MeV, where the dominant uncertainty now arises from the experimental moments. Re- 
summation techniques are also shown to reduce renormalization-scale dependence in the relation between &-quark 
MS and pole mass and in the relation between the pole and IS mass. 



1 Introduction 

Comparison of theoretical and experimental moments Mm, defined by 

Mn^ J ^R{s) , Af = l,2,3... (1) 
a (e+e^ hadrons) 47ra^ . 

-K(S) = , Tpt = — — , [2) 

(7pt 6s 

provides a method for determining the MS quark masses [T . This method, combined with recent BES data [2] 
(particularly in the charm threshold region) and O [af) (mass-dependent) perturbative expressions for the moments 
Mtv and for R{s) in the continuum region [5], has resulted in precision determinations of the MS charm and bottom 
quark masses 0]. 

The MS 6-quark mass determined in 4 contains three major sources of uncertainty; 

1. experimental values of resonance and threshold contributions within R{s), 

2. uncertainty in ag (Mz), which enters both the perturbative series for Mn and the QCD continuum contribution 
to experimental moments, 

3. and theoretical uncertainty associated with renormalization scale dependence within the Mn perturbative 
series. 

In the 6-mass estimates of Ref. these three sources of uncertainty play different roles as N varies, but are 
generally comparable in magnitude. Thus as experimental information becomes more precise, the theoretical uncer- 
tainty devolving from renormalization-scale dependence will become increasingly significant, and without theoretical 
improvement, will be the limiting factor in the precision of 6-mass estimates. 
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Techniques for substantially decreasing renormalization scale dependence have been developed and applied to a 
number of MS perturbative processes including semileptonic b decays, light-quark contributions to R{s), and Higgs 
decays j5j. These techniques use the appropriate rcnormalization-group (RG) equation for each process to determine 
and resum all logarithmic contributions in the perturbation series that are explicitly determined by the RG equation. 
In this paper we extend and apply such techniques to the perturbative series for the 6-quark contributions to Mn, 
effectively eliminating their renormalization-scale dependence as a source of theoretical uncertainty. It should be 
noted that Ref. 0] does not attempt to estimate the uncertainty associated with such higher-order contributions, so 
the uncertainties quoted in Ref. 0] may be underestimated. 

In Section 13 we develop an analysis of the 0(25 MeV) residual renormalization-scale dependence characterizing 
the extraction of the MS &-quark mass from the first four moments of R{s). In Sectional we demonstrate how this 
scale dependence is essentially eliminated upon incorporating the closed-form summation of leading and successively- 
subleading logarithms within the perturbative series for N = {1,2,3,4} moments of 6-quark contributions to R{s). 
This procedure is also shown to reduce theoretical uncertainties associated with the choice of N, as well as leading 
to a modest (14 MeV) elevation in the central value for mb{mb). Finally, by assuming power-law growth in the 
(RG-undetermined) perturbative contributions, the effect of the (unknown) next-order perturbative contributions is 
estimated. These estimated next-order contributions are decreased by the application of resummation techniques, 
providing a further reduction in a source of theoretical uncertainty. 

Renormalization-scale dependence is also shown to exist as a source of uncertainty in the known perturbative 
expressions relating the &-quark MS mass to its corresponding pole mass and IS' mass [3 El- In Section^ we 
explore the scale dependence inherent in the perturbative series relating MS and pole 6-quark masses. In Sectional 
we demonstrate how this uncertainty is resolved by a renormalization-group resummation of this series similar to that 
of Section|31 Finally, in Section|^we discuss the reduction of scale uncertainty via comparable renormalization-group 
resummation of the relationship between the IS and pole 6-quark masses. 



2 Residual Scale Uncertainty of the MS 6-quark Mass 

Following Ref. Kiihn and Steinhauser ^ express the running &-quark MS mass m\j{ii) from moments of the 
&-quark contribution to the experimentally determined electron-positron-annihilation ratio 

Mn'^J^M-s) (3) 
and its (MS) field-theoretical analogue 

Mh'-ijzirr^Y SnM^ (4) 



with the perturbative series Sn related to 11^ (g^) , the 6-quark contribution to the vector-current correlation function, 
via 

If one equates the experimental and theoretical moments, one finds that 

The series Sn is a perturbative series in the QCD couplant x{fi) = as{iJ.)/n, where /i is the renormalization scale 



characterizing M*^ 



N ■ 

oc ] / 2 



3=0 k=a 



where the j = {0, 1, 2} [i.e., up to three-loop] MS coefficients of this series |3] (summarized in Table 6 of Ref.^), are 
tabulated in Table Since Sn has logarithmic dependence on iribip), Eq. @ represents an implicit equation that 
must be solved numerically to determine mb{fJ,). 

^Our coefficients Tj^^ are related to those of Ref. |l|'s Table 6 by Tj^'' = (— l)''c$f ''°'/4. Division by 4 is a consequence of = Q^/i. 
The alternation in sign follows from the argument of our logarithm in Eq. j7J being the inverse of that for the logarithm in the series Cn 
of Ref. E|. 
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0.7898 
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0.4762 


0.8296 


1.1236 



Table 1: Five-flavour {nj = 5) MS series coeSicients for Sn- 



The theoretical expression (@J| for the moments Mn is formally independent of the renormalization scale /i, 

;ity.2 

renormalization-group equation for the series SN^fi): 



as expected for this physically-observable quantity.^ Thus the requirement = dM^j^ /dji^ leads to the following 
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/3o = /?! = Y^, /?2 = 9769/3456 . 



For example, one finds from Eq. JSJ that the coefficients T^^\ ^i^"* ^'^^ '^21'' satisfy the relations 



rW = (/3o-2^)TW 



u ' 

-27V(7i-2)r(^\ 



(8) 

(9) 
(10) 

(11) 
(12) 
(13) 
(14) 



(15) 
(16) 
(17) 



consistent (modulo round-off errors) with the entries in Tabled 

We are interested in exploring both the residual renormalization scale dependence and the iV-dependence of the 
MS benchmark mass nib (mf,) extracted in Ref. since each such dependence is a source of theoretical uncertainty. 
As in Ref. the series ((Tj) for SN{fi) is truncated after its (known) j = 2 terms. Such truncation necessarily 
becomes a source of residual /i dependence. Since we are focusing only on theoretical uncertainties arising from such 
scale dependence, we assume that x{^) four-loop evolves from its Ref. 0] benchmark value x (Mz) = 0.11800/7r to 
a;(10 GeV) — 0.056732, and disregard theoretical uncertainty associated with x {Mz) [and hence a;(10GeV)]. The 
preferred renormalization scale in Ref. 0] for extracting rub {rub) is 10 GeV, and as in Ref. @], the scale dependence 
we consider is over the range 5 GeV < /x < 15 GeV. The extraction of mf, (mf,) occurs first by numerical solution of 
Eq. Q to obtain mf,(/z), and then by evolving mb{fi) downward via Eq. Hll() to the point where n = mf,(/x) = mf, (mf,). 

In Table|21 we list the above-described extractions of nibifJ.) for = 5, 10 and 15 GeV, as obtained from each of 
the A'' = {1,2, 3, 4} moments using values for from Table 7 of Ref. @|. ^ 

The experimental values for Mjy are tabulated in Table 7 of Ref. 
^See Refs. IIUI 1111 and Refs. 1121 1131 for the coefficients in /3(x) and 7(x), respectively. 
"TableElnumbers for = 10 GeV are in agreement with those of Ref. |4|'s Table 8. 
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N = 1 


N ^2 


N = 3 


N = 4 






5 


4.0848 


4.0799 


4.0751 


4.0707 


4.0776 


0.0053 


10 


3.6652 


3.6509 


3.6407 


3.6551 


3.6530 


0.0088 


15 


3.4708 


3.4457 


3.4445 


3.5290 


3.4725 


0.034 



Table 2: (p) as extracted from © via Sn{p) truncated after three-loop terms. All entries are in GeV. The 

bold-face entry [A^ = 2, n= 10 GeV] corresponds to that preferred in Ref. 0] to generate mt {nib). 

The column TOh(/i) of Table El is just the average of rrn,{iJL) taken over the first four moments. The rms spread of 
values over the first four moments is 

1/2 

(18) 

One sees immediately from the table that this rms spread of values for mb(/i) increases dramatically with indicative 
of residual scale dependence. In other words, the error associated with different choices of N is itself a scale dependent 
quantity. 

If Table El represents a valid determination of the quark mass, then for a fixed N , the variation of mh{^) with /i 
should conform with the RG evolution equation Hll|l . However, the /i-dependence of mi,{^) as extracted via Eq. Q 
is not fully consistent with such evolution. In Figure ^ we have plotted the /i-dependence of such extracted values 
for the N — 2 case against the (three loop) RG-evolution following from Eq. (|ll|l for values of /i between 5 GeV 
and 15 GeV To facilitate this comparison, we evolve from the same value TOb(lOGeV) as extracted in Table El (in 
bold). The figure clearly shows a deviation by the /^-dependence extracted from Eq. from that anticipated from 
RG-evolution. Moreover, this deviation becomes progressively pronounced with increased TV. For N = A, Figure 
El shows the plot of extracted versus RG-evolved fi dependence, which exhibits a substantially larger deviation for 
large /i, though a somewhat better fit between fj, — 5 GeV and /i = 10 GeV. The deviation exhibited in the Figures 
of the extracted 6-quark mass from the behaviour expected from RG evolution is significant since it is much larger 
than the effects associated with a change to one-loop higher (or lower) in the RG-evolution curve. 
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Figure 1: Renormalization-scale (/i) dependence of mi,{fj,) extracted via Eq. @ (solid curve) for the N = 2 moment 
compared with the RG evolution of mb{^) (broken curve). For RG evolution, mb(lOGeV) is used as a reference 
scale, and hence the two curves intersect at /x = 10 GeV. 

These identifiable remaining scale dependences, both horizontal (iV-dependence) and vertical (deviation from 
RG-evolution), necessarily percolate into estimates of to^ (m^). In Table Owe display values of rrib {rrib) obtained by 
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Figure 2: Rcnormalization-scale (n) dependence of ^^(/i) extracted via Eq. © (solid curve) for the iV = 4 moment 
compared with the RG evolution of nrn,{^) (broken curve). For RG evolution, mf,(10GeV) is used as a reference 
scale, and hence the two curves intersect sA ^ — 10 GeV. 





N = 1 


N ^2 


AT = 3 


N 




(Tat 


5 


4.2116 


4.2073 


4.2031 


4.1992 


4.2053 


0.0046 


10 


4.2067 


4.1928 


4.1830 


4.1969 


4.1948 


0.0085 


15 


4.2025 


4.1768 


4.1755 


4.2620 


4.2042 


0.035 




0.005 


0.015 


0.014 


0.033 


0.0052 





Table 3: {rrib) as RG-evolved from values listed in Table [3 (i.e., via 3-loop-truncated series Sn)- All 

entries are in GeV. The bold-faced entry [iV = 2,/i = 10] corresponds to preferred choices in Ref. 

evolving Eq. values mf,(/x) extracted from Eq. JHl at the indicated scale fi. The values for = 10 GeV, including 
the preferred N — 2 value, differ inconsequentially (they are 2 MeV larger) from the central values displayed in the 
erratum to Ref. 0, providing a check on our calculation. The column (mf,) of Table |31 is just the average over N 
of mil (mf,) values, as evolved from the indicated choice of /i. 
Similarly, the rms spread over N is just 

1/2 

(19) 

and the scale uncertainty (cr^) is just half the difference between the maximum and minimum value of mt, (mf,) 
over a given column of the table. Thus tJjv is a measure of horizontal (A^-dependence) uncertainty, and is a 
measure of vertical (fJ,—) renormalization-scale uncertainty. For the N — 2 ji ~ 10 GeV preferred case 0, both of 
these uncertainties are indicative of overall 24 MeV theoretical uncertainties in mf, (mf,) that devolve ultimately from 
residual scale dependence in truncating the series Sn- 

Another source of theoretical uncertainty is the effect of next- and higher-order contributions in the perturbative 
series used to determine mf,(/x). Such effects were not considered in the error analysis of 4 . The RG equation is 
capable of determining Tg^^ for n — {1, 2, 3} 

-T^,f = {^^NPl - 2f3,N^ + l7V3^ (20) 

^(2^ = (-8iV^ + 6A^/3o - N(3i -f AN'^^i - 27V/3o7i) ^0,0^ + (-3iV/3o + /3o + 2Af^) t["o^ (21) 
-T^'l^ = (87V + 2iV72 - 87V7i) Tq^^^^ + (27V71 + 2/3o - /3i - 47V) T^'p + (-2/3o + 2N) T^^^^ , (22) 



.N=l 



|m[^'' (mf,) - mf, (mf,)| 



5 



leaving only T-^^^ undetermined. However, we can estimate the effect of these next-order contributions by assuming 

that the approximate power-law growth exhibited by Tq^\ T^'^\ and T^j^^' continues at next-order, resulting in the 
estimates 



(1) 



3,0 



3.64, T, 



(2) 



3,0 



2.15, T, 



(3) 



3,0 



1.42, T, 



(4) 



3,0 



0.784 



In addition to inclusion of the Ts.n, it is necessary to include the next-order [uf = 5) terms 
(33 = 18.8522, 73 = 11.0343 



(23) 
(24) 



in the evolution of the running coupling and mass. These higher-order terms lead to an additional theoretical 
uncertainty of approximately 20 MeV for the N — 2 benchmark m?, (m?,) value, comparable to the renormalization- 
scale uncertainty in Table|3l By comparison, the uncertainty in irib {mb) arising from the experimental inputs (M^^^ 
and as [Mz) ^S]) is approximately 40 MeV, and hence the theoretical and experimental uncertainties are of 
comparable magnitude. 



3 Optimal RG-Improvement of the MS 6-quark Mass 



As noted in the previous section, the higher order Sn{^) series coefficients T^^^"* and are determined via the 



RG-equation (|SJ) from the leading series coefficient Tq^-* . Similarly, the three-loop coefficient ^ can be obtained via 



(N) 



Eq. (|17l) from the one and two-loop series terms T^q' and T^q ■ In fact, the RG-equation |(SJ) is much more powerful 
than any use we have made of it so far. Given the calculated values of T'g^'', T^^^^ and Tj^', one can determine 



respectively every leading-logarithm (LL) coefhcient Tj^\ every next-to-leading- logarithm (NLL) coefficient T'^ jZi, 

and every next-to-next-to-leading logarithm (NNLL) coefficient T^^}^2 the series expansion Q for S'7v(/^)- 

The procedure of optimal RG improvement |S] involves the summation to all orders of leading and progressively 
subleading logarithms within a series, a process that has been seen to reduce significantly the renormalization scale 
dependence in a wide variety of processes 0. For the case at hand, we wish to include every RG-accessible coefficient 
Tj^^ in the series S'Ar(/i) in order to extract via Eq. Q an MS 5-quark mass rrib (mb) that is free (or nearly so) of 
the residual scale dependence evident in Table 13 To do this, we first organize the series ((JJ as follows: 



(N) 



n=0 



(25) 



where 



00 

k—n 



(26) 



with u = xL amounts to an LL summation when n = 0, an NLL summation when n = 1, and an NNLL summation 
when n = 2. If we substitute Eq. (|25|l into Eq. 0, we generate a succession of first-order differential equations for 
these summations: 



2N 



5r(.)=o, 



d 



du 
(1 - Pou 

if3iu - 2 



(2iV - po) 



'{u) 



d_ 
du 
d_ 
du 



{2N - 2/3o 
(A ~ 27V71) 



{Piu - 2) 



du 



2Nji 



(/32"-27i) 



Eq. provides the initial conditions Sq^\o) — tI^^\ 



du 

rpiN) 

1,0 



2Nj2 



AN) 



(0) 



— T2^^ , the set of all known Tj, p 



(27) 
(28) 

(29) 

(AT) 



coefficients of the series S^ifJ^) [see Table^. With these initial conditions, Eqs. ||77I) - (P5|) can be successively solved, 
and the optimally RG-improved series is found to be 



S^[xL] = S^o^\xL) + xS[^\xL) + x'sf\xL). 



(30) 
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For 



w=l~ PoxL, A = -2N/f3o, 
we find from Eqs. (HTj) and that 



and 



S[^^ (xL) = Bw-"" + r| ' -B + C\ogw 



-A-l 



with 



Substituting Eqs. ^ and into Eq. ^ we obtain the solution 



D 



w 



(N) 



D 



H 



with 



D= (3iAB + (32AT^^o^ - iPi - 2N-fi) B + 2Nj2T^"o' /Po 



2 ---£; + F + GlogH + ylog"w; 



(N) 



E : 



F : 
G 

H 



(2/3o - (3i) AB + Pi [(1 + A) (r[^' -B^-C 



+ (27i/?o - P2) AT^^f + [b- t(^)) (/3i - 27V71) 

C(^/3i + 27V7i)//3o, 
(1 + A) {t['1^ -B)-c\ (2/3o - //3o, 
(2/5o-/3i)(l + A)C//3o. 



//3o 



(31) 
(32) 
(33) 
(34) 

(35) 

(36) 

(37) 

(38) 
(39) 
(40) 



The extraction of (m^) now proceeds analogously to that in the previous section, except that the series S'Ar(/i) 
is now in the optimally RG-improved form (|30|l . rather than the truncation 



S 



j=0 k=0 



(41) 



of Eq. O to one-, two-, and three-loop contributions utilized in Section 2 (and employed in Rcf. |4j). For a given 
choice of renormalization scale /i, the values of mb(/i) we extract from Eqs. 10} and H3U|) are tabulated in Table 0] 
As before, the average mb(/x) is over N = {1,2,3,4} values of mb{fj,), and the rms spread over these values is given 
by Eq. p8l) . These spreads are seen to be significantly less than those of Table |21 We thus see that the extracted 
values for mb{fJ-) from different values of N are in much better agreement when Sn is RG-improved. Moreover, the 
substantial increase of ctat with /i characterizing Table|2(and indicative of residual scale dependence) does not occur 
in Table0] In TableQ] ctn is essentially static at 4-5 MeV for fj, between 5 and 15 GeV, corresponding to a small fixed 
theoretical uncertainty associated with the choice for N. Thus, RG-improvement is seen to disentangle (vertical) 
scale-uncertainties from (horizontal) iV-uncertainties, as well as to reduce the magnitude of such A''-uncertainties. 

In contrast to Figs. Q] and [3 the Table Rvalues of mf,(/z) extracted from Eq. 10 via the resummed of Eq. 
(|30|l are fully consistent with the RG-evolution equation Hll() . For N — 2 Figure |2| plots both extracted values of 
mb{n) for /X between 5 GeV and 15 GeV, as well as the values evolved [via Eq. Hll|l ] directly from the extracted value 
to), (10 GeV) = 3.6653 GeV. The points coincide to within the visual resolution of the figure, indicative of purely 
RG scale-dependence for values of rribin) extracted via Eq. H30(l .^ 

^The diflference between extracted and evolved values at /x = 15 GeV is less than 10"* GeV, based on identical starting values at 
/i = 10 GeV. This difference is characteristic of the effect generated by including next-order terms in the RG-evolution equation. 
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N = 1 


N ^2 


TV = 3 


N = A 




(Tat 


5 


4.0851 


4.0797 


4.0750 


4.0717 


4.0779 


0.0050 


10 


3.6698 


3.6653 


3.6610 


3.6582 


3.6636 


0.0044 


15 


3.4793 


3.4751 


3.4712 


3.4687 


3.4736 


0.0040 



Table 4: (/i) as extracted from optimally RG- improved Sn{^). All entries are in GeV. The average mf,(/i) and 
the rms spread over the four values of N are calculated as in Table |2 All entries are in GeV. 



4.1 




3.4 -\ , , 1 1 

5 7 9 11 13 15 

H (GeV) 

Fi gure 3: Renormalization-scale (/i) dependence of mi,(^^) extracted via substitution of the resummed quantity 
into © (solid curve) for the N ^ 2 moment compared with the RG evolution of mi,{fi) (broken curve). For RG 
evolution, mf,(10 GeV) is used as a reference scale, and hence the two curves intersect at /i = 10 GeV. The two curves 
overlap completely within the resolution of the figure. Similar overlap occurs for all other moments considered {i.e. 
iV = 1,3,4). 

The absence of any additional residual scale dependence naturally carries over to the RG-evolution of extracted 
mb(/i) to TTifc (mfc). In Tablc[Slwe list values of nib {nib) obtained by evolution of the value nib{n) extracted at the 
scale n for the indicated moments N. RG-improvement is seen in Table0]to virtually eliminate the scale uncertainty 
(cr^) evident in Table |3| The only theoretical uncertainty still evident is horizontal, the crjv associated with different 
choices of N, and this uncertainty is both small (0(4 MeV)) and static as varies from 5 to 15 GeV. For Ref. @]'s 
phenomenologically motivated choice ^ = 10 GeV, ctat for the RG improved case (TableEl is less than half the value 
for (Tjv when Sn is truncated (Table OJ. 

This virtual elimination of residual scale dependence, as evident in Tabled also leads to a changed central value 
for nib {nib) relative to that of the erratum to Ref. 0]. The central value quoted in the erratum is 4.191 GeV, based 
on choices /i = 10 GeV, N = 2. The corresponding value in our Tabled is 4.193 GeV, and the 2 MeV discrepancy 
is insignificant compared to the erratum estimate of theoretical uncertainty (±51 MeV), or relative to Table |21s 
±15 MeV and ±9 MeV vertical and horizontal theoretical uncertainties associated with the choice of /i and N. The 
corresponding N = 2, = 10 GeV value in Table [S] is 4.207 GeV with aggregate 4.5 MeV horizontal and vertical 
uncertainties. 

Thus, the incorporation of an optimally RG-improved perturbative series SN{fi) is seen to eliminate the renor- 
malization scale theoretical uncertainty and halve the moment-dependence theoretical uncertainty, with a 14 MeV 
increase in the corresponding central value 4.207 GeV. Note the near equivalence of this value with the nib {nib) = 
4.205 GeV values found from the nib (nib) averages over the first four moments (Table EJ. 

As in the previous section, the approximated values H23|) for Tg^-* can be used to estimate the effect of higher-order 
perturbative corrections on the extraction of nib {"mb)- However, with the input of T^^\ it is necessary to extend 
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N = 1 


N = 2 


iV = 3 


iV = 4 




aN 


5 


4.2118 


4.2071 


4.2030 


4.2002 


4.2055 


0.0044 


10 


4.2112 


4.2068 


4.2027 


4.2000 


4.2052 


0.0042 


15 


4.2111 


4.2069 


4.2028 


4.2003 


4.2053 


0.0041 




0.0007 


0.0003 


0.0003 


0.0003 


0.0003 





Table 5: to[ (to&) as RG-evolved from to[ (/i) values listed in Table 0] {i.e., via the RG-improved series S^). 
rnfc (mfc), the average of rrih [nii,) over N, aN, the rms spread over N, and the scale uncertainty ct^ are obtained as 
m Table El All entries are in GeV. 



(|30|l to include an slf^\xL) term 



S^[xL] = S^„^\xL) + xS[^\xL) + x^si^\xL) + x^si^\xL) 



(42) 



The expression for S^^\xL) is determined by an extension of equations H27|l -H29 |l . and the final result can be 
extracted by appropriate modifications of .5, . The resulting expression is 



AN) 



(xL) 



K 



-w 



Rw 



-A~2 



M N 

T ~ T 

log'(u)) + 



w 



-A-1 



N 



-w 



-A-l 



\og{w) + {P-Q + 2R)w 



-A-2 



{Q ~2R)w-'^~^\og{w) 



K M N 



Uw 



-A-3 



3 2 4 

.2. ^ Y 



P + Q-2R + T3^o]w 



log(w) + —W log (w) + —W 



2 ° ' ' 3 

where the (A^-dependent) coefficients are given by 



K 

M 



A 



2iV73ro^,o^ + (-27V72 + P2)B + {~2Nji + 2/5i) 



D 



272 



/3c 



(N) 



0,0 



271 - 



(-2A^72 + f32)[B- 0' ) + (-2A^7i + 2/3i) {F - E) 



B 



270 - 



l3o 



A 



1 

Jo 



(t^^o^ -b'^{1 + A)-c]^ + [Eil + A)- F{2 + A)] ^ 



/3o 



N = {(A/32 + 2iV72) C+[{A- + 2N-fi] F} — 

Po 



p = 27; 



f32 



/3o 

(-2A^7i + 2f3i 



Pi 



2^0- f^] [{I + A)E~ {2 + A)F] 



D 



' Po 



G-i2 + A)[T^^,^-^-E + F 



(l + A)- [(-2A^7i + 2/3i) G+iH-{2 + A)G)Pi] — 

Po 



i?=(/3iA + 2iV7i)-^ 
^Po 
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270-^ 

Po 



(2 + A)(r(^) 



u 

V 

>^-f27o-|)(2 + A)^ 



2jo~^][{2 + A)G-H] 
Po, 



H 



(43) 



(44) 
(45) 

(46) 

(47) 

(48) 
(49) 
(50) 
(51) 
(52) 



The effect of these higher-order terms lead to an approximate 7 MeV uncertainty in the N = 2 benchmark m?, (mf,) 
value. This is a significant reduction compared with the 20 MeV uncertainty occurring for the un-summed case 
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described previously. By comparison, the uncertainty in the resummed mt {rrib) arising from the experimental inputs 
{M2^^ |3j and as{Mz) Ql]) is approximately 40MeV, and hence the resummation analysis reduces theoretical 
uncertainties to a level well below the experimental uncertainties. 

4 Residual Scale Dependence of the Pole Mass 

The series T(/i) — T [x{fi), L{fi)] relating the (RG-invariant) pole &-quark mass and the MS mass mb{n) is given 
by 

MP ^mbi^I)T[xifi),Li^i)], (53) 



r[x,L] = l + ^^T,,fea;^L\ (54) 

j=0 k=0 

where x{iJ,) = as(/i)/7r and L{ii) — log {^^"^ / m1{iij) as before, and where the known series coefficients in the _L = 
limit are |H] 

Ti,o = 4/3, 

r2,o = -1.0414 (n/ - 1) + 13.4434 9.2778 

Tg^o = 0.6527 (n/ - 1)^ - 26.655 {uf - 1) + 190.595 94.4182 (55) 

If one wishes to relate the pole mass to nih (mf,) (i.e., to the point /i = ■mb{p)), all logarithms in the series 1)54(1 are 
zero and only the coefficients Tj Q contribute. 

Since is a RG-invariant quantity, renormalization scale dependence can be studied by explicitly varying 
rubip) through the RG equation as /i is varied. Given the residual scale dependence implicit in any truncation of 
the series (|54|l . one can argue that the choice of scale fi, generally motivated by experimental information [as is used 
to determine M'^^ via Eq. Q], should be consistently maintained. For example, the benchmark mass mf, (mf,) is 
obtained in Ref. ^ via a determination of TOf,(10 GeV), and then by subsequent evolution [Eq. I|ll() ] to the point 
/i = mf,(/x). 

A measure of the residual scale-dependence implicit in the determination of the pole mass from the MS mass 
obtained via Ref. 4, methodology would be the difference between 

1. the /i = 10 GeV pole mass obtained via Eqs. H53|) and (|54|l by incorporating within the logarithm L{fi) the 
value TO;,(10 GeV) actually extracted from Eq. ©, and 

2. the ^ = mi, (mi,) pole mass obtained via Eqs. (|53() and (|54|l from similar incorporation of mi, (mi,), as evolved 
via Eq. Hll|l from the extracted value mf,(10 GeV). 

We emphasize that if the input values mi, (mi,) and mb(10 GeV) are RG-consistent, then any discrepancy between 
obtained by these two procedures is a reflection of residual scale dependence arising solely from the truncation 
of the series (|54|l . 

To implement this comparison, we need to know the coefficients Tj,k with k ^ for j — {1,2,3}. Since the 
pole mass is an RG-invariant, dM^ /dfi^ = 0. One then finds from Eq. (|53|l the following RG-equation for the series 
r[x,L]: 



(l-27(x))^+^(^)^+^(^) 



T[x, L] = (56) 



If one substitutes the series ((54|l into the above equation, and then utilizes the series pi(l and H13|l for 7(2:) and (3{x), 
one finds after a little algebra that 

Tm = 1, T2.2-(l + /?o)/2 (57) 
T2,i = (1 + l3o) Ti.o + 71-2 5.4208 (58) 

^ , , , rif — 5 1015 , 

^3.3 = (1 + /3o) (1 + 2/3o) /6 ^ (59) 
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T3,2 = (1 + f3o) (1 + 2/3o) Ti,o/2 + /?o (71 - 3) + 71 + ^ - 2 

2 (60) 

13.1053, 

T3,i - (1 + 2/3o) T2,o + + 71 - 2/3o - 2) ri,o + 72-471+4 

(61) 

42.3366 

We then can employ the three-loop series 

3 3 

r[x,L]=^^r,,fc:r-'i'= (62) 

j = l k=0 

with series coefficients given in Eqs. (|55|) and l|57|l - H61|l to compare the pole mass obtained from the extracted 
value m6(10GeV) to that from the correspondingly RG-evolved value irib (mb). To be consistent with having three 
subleading orders in x in series (jSU, we utilize Eqs. (fTT|l and l(T^ to four-loop order to evolve x(/x) and mf,(/z) 
[/33 = 18.8522 11 , 73 11.0343 13 ] from the same a:(10 GeV) — 0.056732 reference couplant value [as evolved from 
an assumed a, {Mz) = 0.11800] used throughout. Using the ^l = 10 GeV, N = 2 value mf^(lOGeV) = 3.651 GeV 
of Table 12 as a springboard value, and its corresponding value mi, (mi,) = 4.19 GeV (Table 0), we find somewhat 
different pole masses for different choices of /j,: 

/i = 10 GeV, AfP = 4.82 GeV, (63) 
/i = 4.19 GeV, MP = 4.94 GeV, (64) 

indicative of 120 MeV residual scale uncertainty. We emphasize that this uncertainty arises entirely from the trun- 
cation of the series (|54f) . and is independent of scale uncertainties in the extraction of r7if,(/i). Had we used the 
corresponding N — 2, ^ — 10 GeV values mf,(10GeV) = 3.665 GeV, mb{mb) — 4.21 GeV obtained in Tables 4 
and 5 via optimal RG improvement, the corresponding pole masses still exhibit virtually the same residual scale 
uncertainty: 

/i = 10 GeV, MP = 4.84 GeV, (65) 
fi = 4.21 GeV, MP = 4.95 GeV. (66) 

Consequently, there appears to be a surprisingly large 110-120 MeV theoretical uncertainty implicit in the determi- 
nation of the pole mass from the MS mass mb (fj.) , an uncertainty devolving ultimately from truncation of the series 
(I54|l after its O (x^) terms. In the section which follows, we will optimally RG-improve the series (|54ll to reduce this 
residual scale uncertainty by a factor of 15. 



5 Optimal RG Improvement of the Pole Mass 

Optimal RG- improvement of the series T[x,L] follows along the same lines as described in Section 3 for the series 
Sn[x,L]. We express the series (|54|) in the form 

00 

T[x,L] = Y,Tn{xL)x" (67) 

where 

?^«(u) = ^Tfc,fc_„u'=-" (68) 

k—n 



with Tq{xL) encompassing the LL summation (all values of T^^k), Ti{xL) encompassing the NLL summation, etc. 
Since T3,o is known [Eq. ^^], the functions Tq(xL), Ti{xL), T2{xL) and the N^LL summation T^{xL) within Eq. (|^ 
are all RG accessible. If we substitute Eq. H67I) into the RG-equation H56() . we find the following set of sequentially 
solvable first order differential equations for Tn{u), n = {0,1,2,3}, with initial conditions T'„(0) = Tnfi [Eq. H68|l ] 
explicitly listed in Eq. (|55|l : 



(1 - M — 
du 



Toiu) = 0, (69) 
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(1 - M ^ - (1 + /3o) 



Ti{u) 







(1 - /3o«) :r - (1 + 2/3o) 



(/3i«-2) — + (7i+/3i) 



ri(u) 



= - 

+ 



(/Jm - 2) 1^ + (71 + 2/3i) 



T2(u) + 
To(u), 

TsW + 

Ti{u) + 

Given u = xL and the definition w = 1 — poxL, the solutions to the above four equations are 
To{xL) = w-^, 

Ti{xL) = Bvj-^ + w-^-^ [Ti,o - B + Clog(w)] , 



(/32M-27i) — +72 
au 

(1 - pou) 1^ - (1 + 3/3o) 



(/32W-27i) — + (72 +/32) 



(/?3U - 272) ^ + 73 



T2(U) 

Mu). 



-A-2 



-^-1 [E - F + F\og{w)] 
D 



T2fl---E + F + G\og{w) + y log2(«;) 



TsixL) =yW-^ + w-^-' \^—-- + - log(«;) j 

+ [P-Q + 2R+{Q- 2R) \og{w) +R\og^{ 



-A--i, 



K M N 

■— - — + — -P + g-2i? + T3,o + C/logH 



V Y 
— log^{w) + — log^('u;) }, 



where parameters A -Y are now given by 

A = l/po, B = (A/3i-7i)//3o, C = (2/3o - /3i) A//3o, 
D = [piAB + P2A - + 71) S - 72] //3o, 
/3o^ = (2/3o - + [(1 + A) (Ti,o - - C] /3i 

+ (2/3o7i - /32) ^ + (B - ri,o) (A + 71) , 

F=(A/3i-7i)C//?o, 

G = [(1 + A) (Ti.o - i?) - C] (2/3o - //3o, 
if = (2/3o-/3i)(l + A)C//3o, 

/3oi^ = [A {133 + P2B + (3iD/2) - B (72 + /32) - 73 + (71 + 2/3i) £'/2] , 
0oM = [(2/3o72 - /33) + B (2/3o7i - P2) + D (2/3o - /2] A 

+ [72 + /32 - (A + l)/32] {B - Ti,o) - C/32 

+ (71 + 2/3i) (F - £;) + [E{1 + A)- F{2 + ^)] 

TV = [C (A/32 - 72) + [(A - l)/3i - 71]] //3o, 

/3oP = (2/3o7i - /32) [(1 + A) (Ti,o - B) - C] 
+ (2/3o - /3i) [(1 + A)E - (2 + 

- (71 + 2/3i) (T2,o - £»/2 - £ + F) 

- A [G - (2 + A)] (T2,o - r'/2 - F + F) , 
/3oQ = [(2/3o7i - /32) G + (2/3o - Pi) F] (1 + A) 

-(7i + 2/3i)G-[iJ-(2 + A)G]/?i, 
i?= (AA-7i)i//2/3o, 

t/ = (2/3o - [(2 + A) (T2,o - D/2 - F + F) - G] //3o, 



(70) 
(71) 

(72) 



(73) 
(74) 

(75) 



(76) 



(77) 
(78) 

(79) 

(80) 
(81) 
(82) 
(83) 

(84) 

(85) 

(86) 

(87) 

(88) 
(89) 
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V =[2(3o-Pi)[{2 + A)G-H] //3o, 
y = (2/3o-/5i)(2 + A)F/2/3o. 



(90) 
(91) 



We now can compare the pole mass obtained via Eq. H53|) from the first fom- terms of the series (|67|l . which 
include all RG-accessible coefficients, 

T[x, L] = ToixL) + xTi{xL) + x^T2{xL) + x^T^ixL), (92) 

to the pole mass obtained via the three-loop series Ht)2|) . If we utilize Table Els N = 2, ^ — 10 GeV extracted value 
TO6(10 GeV) = 3.665 GeV and incorporate Eq. ^ into Eq. we find that 

fi = 10 GeV, MP = 4.959 GeV (93) 
fi = 4.207 GeV, = 4.951 GeV (94) 

This 0{8 MeV) uncertainty is a remarkable improvement over the 0(110 MeV) uncertainty [Eqs. (|65|l and H66|l ] 
following from the same input assumptions, but with T[x,L] given by the three loop series (|62|1 . Thus, the RG- 
improved series H92|l removes virtually all the residual scale dependence in the relation between the pole and MS mass. 
This is corroborated in Fig.01 in which the pole masses obtained via the three loop series (|62|l and the RG-improved 
series (|92|) are compared directly, given identical m;,(10GeV) = 3.6636 GeV anchoring values for the MS mass [the 
/i = 10 GeV value for rrTb{rnb) in Table The pole mass for the RG-improved case exhibits very little dependence 
on /i. Thus optimal RG-improvement of the perturbative series (|54|l is seen to remove virtually all of the substantial 
(0(110-120 MeV)) residual uncertainty characterizing the relationship between the MS and the pole 6-quark mass. 




Figure 4: Renormalization scale dependence of the resummed expression of the pole mass (broken curve) compared 
with the unsummed expression (solid curve). The MS mass mb{fi) used as input to this comparison is based on RG 
evolution from the reference value m(,(10 GeV) = 3.66 GeV as outlined in the text. 
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Optimal RG Improvement of M^^ /Mj^ 



The 6-quark 15 mass (^M^^), defined to be half the perturbative mass of a (theoretical) 35*1 bb meson, has been 
determined via relations between masses and widths of T mesons to moments of the 6-quark vector current correlation 
function 14 . This 15 mass is related to the pole mass Mj^ (which is much more sensitive to Aqcd) via the 
perturbative relationship |7||S|, 



Ml 



IS 



Ml 



1 



2n^x^{^i) 



for which the series W{^) is known in full to two subleading orders: 

3/i 



w{^Ji)^w[x{^i), i[iJi)] , %) = iog 



ATTx{fl)MP 



(95) 



(96) 
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n—l m— 



where 



•^1,0 = -Tr 



97 lln/ nf=5 181 
6" 



9 



18 



T2.0 = 33 7.95 - 40.96571/ + 1.1629n^ 162.19 

n/— 5 23 



rif — 5 



T2,l =6/3oTi,o-8/?2+4/?] 



n/=5 1151 
* 12 



T2,2 = 12/3o' 



2 nf—d 529 

^ T2"' 



(97) 

(98) 
(99) 

(100) 

(101) 
(102) 



In principle, M^^ and are both RG-invariant entities independent of the renormahzation scale /i. Consequently, 
one can obtain the following RG-equation for VF[a;,^] by requiring that [x^(/x)W^(/i)] = 0: 

w[x,e] = o, (103) 



2/3(x)+x/3(x)— + [|-/3(a;) 



d_ 

d£ 



with (3{x) given by Eq. (|13|l . If one substitutes the series (|97|l into the RG-equation p03|) . one easily corroborates 
the results ^U^-^U^. 

Optimal RG improvement of the series is obtained by expressing the series in the form 



W[X,£] = J2^"^n{xi) 



n=0 



where 



Wnixi) = ^rfe,fc^„(x£)'=- 



(104) 



(105) 



k—n 



As before, the series WQ{xt) is inclusive of all LL coefficients Tk^k in the series (j97|l . Similarly, Wi{x£) includes all 
NLL coefficients r^.k-i, and W2{x£) includes all NNLL coefhcients Tk,k-2- Upon substituting Eq. I|104|) into the RG 
equation 1103|l . we obtain successive first order differential equations for the n = {0, 1, 2} cases of series (I105|l : 



(1 - 2/3ow) — - 4/5o 
au 

(1 - 2(3ou) - 6/?o 
au 



Wo{u) = 0, 



Wi{u) 



2 {(3iu -!3q) — + 4pi 
au 



Woiu), 



(1 - 2/3o) 



du 



W2iu) = [2 {piu - l3o) + 6/3i] Wi{u) 

+ [2 {(32U-Pi)+AP2]Wo{u). 



(106) 
(107) 

(108) 



Given the initial conditions VKo(O) = tq^ = 1, Wi(0) = ti.q and 14^2(0) = T2_o evident from the definition (|105|) of 
W„(m), we obtain the following solutions to Eqs. H106|l . (|107|l and H108|l : 



Wo{xe) = (1 - 2/3ox£y 
Wi{x£) = (1 - 2(3ox£y 



Tuo + 4 /3o 



2/3o 



log (1 - 2f3ox£) 



W2{x£) {1 - 2PQx£y^ 

[t2,o -U + VX log (1 - 2I3qx£) + SV^ log^ (1 - 2/3oa;£)] 



(1 + 2l3ox£y 



where 



U^2 



'01 _^ 
A. 



(109) 
(110) 

(111) 
(112) 
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v = 

X 



(A - 2/32) //3o, 
-3ri,o + 4/3o-2/3i//3o, 



(113) 
(114) 



with Ti^o and r2,o as given in Eq. pU2|) . 

Let us first consider any residual scale dependence in the relation (|95|l arising from truncation of the series (|97|l 
after the known coefficients (I1U2|I . It has been argued in Ref. [S] that the relation (|95l) is to be utilized at soft 
momentum scales (1.5 GeV < /i < 3.5 GeV), because the 15 mass is defined purely from nonrelativistic dynamics. 
Moreover, a soft scale necessarily follows from the non-perturbative nature of the logarithm i\ if ^ is hard, then x 
is small and logarithms are large. Given Eq. H93() 's 6-quark pole mass of 4.96 GeV, for example, we find for the 
truncated series that the 15 &-quark mass varies between 4.58 and 4.73 GeV as /i increases from 1.5 and 3.5 GeV, a 
150 MeV residual scale uncertainty. This range is substantially diminished if we utilize the RG-improved series 



(115) 



as determined by Eqs. H1U9|) . H11U|) and (|lll|l . for the same region of [i. For the RG improved series, we find that 
the IS* 6-quark mass corresponding to a 4.96 GeV pole mass varies between 4.66 GeV and 4.61 GeV as /i increases 
from 1.5 GeV to 3.5 GeV. Thus, optimal RG improvement reduces the residual scale uncertainty from 150 MeV to 
50 MeV. These results are presented graphically in Fig. [S] 



Ji 935 - 




Figure 5: Renormalization scale dependence of the resummed expression of the ratio M^^ /Mj^ (broken curve) 
compared with the unsummed expression (solid curve) . The range considered for /i is the "soft" region advocated in 
[H], and = 4.96 GeV is used as an input value. 

It is, of course, more realistic to extract a 6-quark pole mass from a phenomenological determination of the IS* 
6-quark mass. Using the central value 4.71 GeV for the IS mass "9 , one can invert the relation H95|l numerically for 
both the truncated and the RG-improved versions of the series (19 7|) . The results of this inversion are displayed in 
Fig. El and are indicative of a pole mass somewhat above 5 GeV. Once again, however, an 0(150 MeV) theoretical 
scale uncertainty for the truncated case is reduced to a 40 MeV scale uncertainty using the RG improved series. Note 
that the crossing point in both figures is the soft-fi point at which the logarithm £ is equal to zero. It is evident from 
the initial conditions for Wq, Wi and W2 that the RG summed series and the truncated series are equivalent at this 
point. 



7 Conclusions 

As demonstrated in Section[21 the procedure for extracting the MS 6-quark mass from empirical moments Mjv of R{s) 
necessarily exhibits theoretical dependencies on the choice of renormalization scale (fj.) , the choice of moment (N) , 
and the effect of higher-order perturbative contributions. Omitting coupling-constant and experimental uncertainties, 
we have shown that an analysis based upon the Ref. U choices /i = 10 GeV and N = 2 leads to the extraction of an 
MS 6-quark mass 

mt (mfc) = 4.193 GeV ± 15 MeV ± 8.5 MeV ± 20 MeV. (116) 
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Figure 6: Renormalization scale dependence of the resummed expression for (broken curve) compared with the 
unsummed expression (solid curve). The input value M^^ = 4.71 GeV and the range considered for /x follows from 
Ref. 0. 



The first theoretical uncertainty is associated with a ±5 GeV variation of the renormalization scale n @], the second 
reflects the moment dependence in letting the choice of moment N vary from 1 to 4, and the third is an estimate of 
higher-order perturbative contributions. 

In Sectional the perturbative series from which this prediction is obtained is optimally RG- improved via the all- 
orders summation of that series' leading, next-to-leading, and next-to-next-to- leading logarithms, i.e., the summation 
of all RG-accessible logarithms in the perturbative series. In addition, an estimated value for T^^^ based on the 
approximate power-law growth of the RG-undetermined perturbative coefficients allows an all-orders summation to a 
further subleading-logarithm order. If input assumptions leading to Eq. (|116|l are otherwise unchanged, the optimally 
RG-improved MS mass is then found to be 

mb (TOfc) = 4.207 GeV ± 0.3 MeV ± 4.2 MeV ± 7 MeV. (117) 

As before, the renormalization-scale and moment-dependence uncertainties displayed above are associated respec- 
tively with varying fi by ±5 GeV, with varying N from 1 to 4, and with the effect of higher-order perturbative 
contributions as estimated via (|23|) . 

The reduction in these latter uncertainties associated with the choice of moment N and the (estimated) higher- 
order perturbative contributions is an unanticipated but welcome feature of the RG-summation developed in Section 
13 Indeed, prior to such improvement, the uncertainty (ittv) devolving from varying N is seen in Table to increase 
quite drastically with /z (ctat = 35 MeV at /i = 15 GeV). After RG- improvement, however, this iV-uncertainty is 
reduced to 4 MeV levels regardless of the choice for fi (TablclSJ. Since the resummation analysis reduces the theoretical 
errors to levels negligible compared with those of experimental inputs {M^^ P' and as (Mz) = 0.119 ± 0.002 |15|). 
the benchmark mi, (mt) determination 0] is modified to 

TOfc(TOfc) = 4.207 GeV ±40 MeV . (118) 

With the reduced theoretical uncertainties and with use of the range as (Mz) = 0.119 ± 0.002 ^H], the dominant 
source of uncertainty now arises from the experimental moments M'^^ . 

Renormalization scale dependence inherent in relations between the 5-quark pole mass and corresponding b- 
quark MS and IS* masses is shown in Sections |5l and to be similarly reduced via optimal RG-improvement of 
the perturbative series characterizing such relations. Comparison of "RG- unimproved" extractions of the pole mass 
[Eqs. 1)65(1 and 1)66(1 ] to RG-improved extractions [Eqs. 1(93(1 and 1(94(1 ] for which all input information is otherwise 
equivalent indicates a reduction from 110 MeV to 8 MeV in the variation of the pole mass with renormalization scale 
as that scale varies from 4.2 GeV to 10 GeV. Moreover, for the improved case, the central-value pole mass is found 
to be near the high end of the range for the "unimproved" pole mass. A similar elevation with RG-improvement 
characterizes the pole mass extracted from the IS* mass, as discussed in Sectional In Fig.El the summation of RG- 
accessible logarithms is shown to lead to a less scale-dependent and somewhat larger pole mass extracted from an 
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assumed 4.71 GcV 6-quark IS-mass than would occur in the absence of such RG-improvement. In particular, there is 
a reduction from 140 MeV to 40 MeV in the variation of the pole mass as the renormalization scale /z varies between 
1.5 GeV and 3.5 GeV. To our knowledge, such renormalization-scale theoretical uncertainties have not previously 

been considered in detail. 

The important point common to all of the cases considered above, however, is that (often-ignored) theoretical 
uncertainties necessarily follow from the residual renormalization-scale dependence characterizing the truncation of 
phcnomcnological pcrturbativc scries, and that such uncertainties may be substantially reduced, if not eliminated, 
by improving such series to include summation of all higher order RG-accessible contributions. Such resummation 
techniques should thus prove to be of increasing value as phenomenological and experimental inputs into 6-mass 
determinations become more precise. 
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